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This talk summarizes the results of recent calculations of how fluctuations within the 
solar medium can influence resonant neutrino oscillations within the sun. Although 
initial calculations pointed to helioseismic waves as possibly producing detectable 
effects, recent more careful calculations show this not to be true. Those features of 
C***~ ' fluctuations which maximize their influence on neutrino propagation are identified, 

Q\ ' and are likely to have implications for supernovae and the early universe. 

> : 

^1 Introduction 

O-JDhe results from solar-neutrino experiments JIJ now strongly suggest that the experimentally-observed 
"^olar- neutrino spectrum differs from the predictions of astrophysical theory |2| . The observed spectrum 



^Qs| consistent with a virtually complete elimination of the flux of intermediate-energy 7 Be neutrinos, 
together with a partial suppression of the higher-energy 8 B neutrino flux. This kind of pattern is 



^--notoriously difficult to obtain using astrophysical variations on the solar-model theme (notwithstanding 
^the recent proposal of significant 3 He convection ||), suggesting instead an interpretation based on 
Nonstandard neutrino properties. Perhaps the most successful such neutrino-based explanation is the 
QvISW mechanism [|j], which proposes the resonant conversion of neutrino flavours within the solar 
>medium, providing an excellent description of the experiments using theoretically-plausible neutrino 

•rH , 

KN)arameters. 
H 

The emergence of the MSW mechanism as the leading explanation for the solar-neutrino problem 
has motivated more detailed scrutiny of the approximations on which its predictions rest. For example, 
the original MSW analysis starts from a mean-field treatment of the background through which the 
neutrinos propagate. More recent work has since investigated how corrections to the mean-field picture 
might influence the neutrino survival probability || H, S H H3- These studies indicate that there 
are astrophysical situations for which fluctuations might significantly influence neutrino propagation 
despite the extreme weakness of neutrino interactions. The two different ways in which they might do 
so are: 

1. Fluctuations can scatter neutrinos away from the forward direction, and so shorten the effective 
neutrino mean free path a the medium. 

2. For situations where the mean particle densities cause resonant oscillations amongst neutrino 
flavours, fluctuations generically act to ruin the quality of the resonance. 



Explorations of the practical implications of these effects are in their infancy. Naturally enough, 
the first applications have been to solar neutrinos, for which the most experimental information is 
available.^] They have therefore focussed on the second of the alternatives listed above, since the 
density of the solar medium precludes seeing the consequences of direct neutrino scattering from 
fluctuations. 

By contrast, MSW oscillations, if they occur, hold out the most promising way of seeing the 
effects of solar fluctuations on neutrino propagation. This is because even though the fluctuation 
effects are small, resonant oscillations make neutrinos more sensitive to fluctuations than they would 
otherwise be. Most interestingly, independent preliminary calculations |7], |8| indicated that reasonbly- 
sized fluctuations in the solar interior might be detectable through the deterioration of the MSW 
resonance which they induce. 

The next sections contain a summary of these estimates as well as subsequent, more precise, 



calculations |I(J of the effects of known solar fluctuations on neutrino propagation. In so doing we follow 
the presentation of refs.|| and |1(J, with which we are most familiar. §2 summarizes which features 
fluctuations must have in order to appreciably affect neutrino oscillations. For solar applications, 
helioseismic waves emerge as the fluctuations with the best chance of spoiling the MSW effect. §3 
then briefly describes the results of detailed simulations of neutrino oscillations in the presence of 
helioseismic waves, culminating with the conclusion that no known source of solar fluctuations is large 
enough to have a detectable influence on solar neutrinos. 



2 General Features 

Much is known about how particles propagate through materials, and general techniques have been 
developed to describe this propagation. This section describes the adaptation of these ideas to neutrino 
physics, as developed in ref. 0. 

2.1 The Origin of Fluctuations 

Since neutrinos interact so weakly with the solar medium, they are negligibly scattered by microscopic 
effects such as thermal fluctuations. As a result, the time evolution of the density matrix, p, describing 
the neutrino state is given by the usual quantum-mechanical expression: 

p(t,n) = U(t,n)p(0)U\t,n), (1) 

where t denotes time, U(t, n) is the evolution operator, and n generically denotes all of those features 
of the solar medium (like the local electron density profile) on which neutrino evolution can depend 
1 See, however, ref. O] for a discussion of some fluctuation effects within supernovae. 



parametrically. For neutrino oscillation experiments we take p to be a matrix in neutrino-flavour space, 
and U is the solution to the evolution equation 

dU 

^ = -iV{t,n)U (2) 

where V(t, n) is the interaction Hamiltonian which couples the neutrinos to each other and to the solar 
medium. For example, for an ultrarelativistic neutrino having three-momentum k and mass matrix m 
propagating through the sun, the dominant flavour- dependent interactions are V = V VSuC + V mat (t,n), 
with: 

V Rj 1- 

' vac ~ „ , i • • • 

2k 

V m&t {t,n) = V2G F g e n e (t), (3) 

where g e = diag(l, 0) is a matrix describing the charged- current coupling to the electron density n e . 
For simplicity we consider here only the case of two active neutrino species. 

For the present purposes fluctuations arise when the electron density, n e , is not constant in space 
or time, since then the solar properties as seen along the trajectories of successive neutrinos can vary. 
The response of a detector to many such neutrinos is then described by an average over the ensemble 
of densities which are seen by the individual particles, weighted by some probability distribution, p(n): 



(0)(t) = J dnp(n) Tr [(D p(t,n)] = Tr [o (p(t))]. (4) 



The physics of the fluctuations is encoded in the probability distribution, p(n), which we define 
in terms of the expansion of n e in a complete set of basis functions, 4>(t), that are assumed to describe 
uncorrelated variables. That is: 



n e (t) = (n e (t)) 



(5) 



where the coefficients Cj, are assumed to be uncorrelated random variables which are Gaussian dis- 
tributed, with vanishing mean: (Cj) = and (CjCk) = o^-b^. 

The most commonly-used choice of basis functions |6|, 0, § are localized in space, corresponding 
to fluctuations which are uncorrelated from point to point. Since this choice is not appropriate for 



helioseismic waves, for these we instead follow refs. || |10| and expand the electron density in terms of 
a complete basis of helioseismic modes [|12|| . 

The brute force approach is to directly solve these equations for U(t,n) and p(t,n), to find the 
electron-neutrino survival probability for any one neutrino after its transit through the sun, and to 
then numerically average the result over the appropriate ensemble. Before describing the results so 
obtained, we first pause to describe an analytic solution which applies in the limit that the correlation 
lengths describing the fluctuations are much smaller than all other scales in the problem. 



2.2 Short Correlation Lengths: The Master Equation 

It pays to recognize that the neutrino-oscillation problem in a random medium necessarily involves 
(at least) two important time scales. This pays because the time-evolution greatly simplifies if one of 
these scales should be much smaller than the other. All of the calculations of solar fluctuations, apart 



from ref. [I(J, are based on this simplifying assumption. This section describes how this approximation 
works. 

One of the basic scales for neutrino oscillations is the oscillation time, r osc , which we define to 
be the time required for different neutrino flavours to acquire significantly different phases as they 
evolve through the medium. By definition, oscillation phenomena always involve evolution over scales 
t > t osc , which cannot be computed perturbatively in the interactions which distinguish between the 
oscillating neutrino species. 

The randomness of the interaction between the neutrino and its environment, described above, 
introduces the other important timescale. More precisely, denoting 5V = V — (V), suppose that the 
correlation (8V(t) 8V(t')) is negligible whenever |t — t'\ is greater than some characteristic scale, r c . 
In this case the neutrino density matrix at time t, p{t), retains a memory of the fluctuations through 
which it has passed only over the time interval (t — r c , t). 

The simplification occurs if r c <C r osc . In this case, the evolution from an initial state, p(t), to a 
state at a later time, p(t + At), can be computed in perturbation theory: 

Ap = p(t + At)-p(t) = ^At, (6) 

so long as At <C t osc . Here Dp/Dt defines a coarse-grained time derivative, which may be computed 
as an explicit convolution of p and V over the time interval (t, t + At). 

On the other hand, if At can also be chosen to satisfy At ^> r c , then Dp/Dt, evaluated at time 
t, depends only on the initial value of p(t), and has no memory of how p has evolved over earlier times. 
This removes the complication of having Dp/Dt involve a convolution over the entire time interval, 
(t, t + At). For example, as applied to the neutrino density matrix in two-by-two flavour-space, the 
twin conditions r c <C At <C r osc lead to: 



Dp 
Dt 



V vac + V MSW (t),p -2G F A(t) (g e Yp + p(g e Y-2g e pg e +0(V 6 ). (7) 



where V MSW (t,n) = \/2G F g e (n e (t)) and the coefficient A(t) of the second term is the correlation 
integral 

A(t) = f dr (5n e (t)5n e (r)) (8) 



which represents fluctuation effects. Notice that the first term on the right-hand-side of eq. (0) describes 
the usual MSW evolution of neutrino flavour. 



Now comes the main point. Eqs. (|6]) and (^) — the Master Equations — define a Markov-like 
process which gives dp(t)/dt ~ Dp/ Dt purely in terms of p(t) and V(t). This relation may be directly 
integrated to obtain p{t') for any later t' . The result holds even for \t' — t\ ^> t osc , even though Dp/ Dt 
is only computed perturbatively in V. 

For two-flavour neutrino oscillations this master equation, including fluctuations, may be solved 
analytically subject to the same assumptions as are typically used to solve the MSW oscillation problem 
|§. The resulting expression for the electron- neutrino survival probability is a generalization of the 



well-known Parke formula [13[ to include the effects of fluctuations: 



P e (t,t') = l+(l-Pj) A<"os2fU/) <-os2tf,„<n. 



(9) 



Here P 



exp 



_tt / sin 2 26V ^ ( $rn 2 h 
2 V cos 26V / V 2fc 



is the 'jump' probability, where h is the scale height for (n e (t)} } 
and 5m 2 is the squared-mass difference between the two neutrino mass eigenstates in vacuum. 9 V is 
the vacuum mixing angle, while 6 m (t) is the matter mixing angle evaluated at the position occupied by 
the neutrino at time t. The contributions of fluctuations are summarized by the coefficient A, which 
is given by 

A = exp [-2 G 2 F f dr A{t) sin 2 20 m (r)l . (10) 
L Jf 

An important consequence of this equation is its implication that the probability of v e survival depends 
almost exclusively on fluctuation properties at the position of the MSW resonance, since sin 2 29 m is 
typically sharply peaked at this point. 

Figure 1 plots the survival probability vs neutrino energy which is predicted for two particular 
models of solar density fluctuations. Successive curves in this plot represent fluctuations having dif- 
ferent values for their correlation length, £, and their relative amplitude, e 2 ~ (5n e 5n e ) / (n e ) 2 at the 
position where the MSW resonance occurs. 

Figure 1 
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Figure 1: The electron-neutrino survival rate vs neutrino energy as predicted by the Master-equation 
method using a particular model of solar fluctuations, and representative solar-neutrino properties. Each 
curve represents a fluctuation having a different value for the one relevant combination of amplitude, e, 
and correlation length, I. See ref. j|| for details. 



2.3 Long Correlation Lengths 

The generalized Parke formula provides an excellent description of neutrino oscillations in media 
having short-correlation-length fluctuations, as has been verified by comparing it with direct numerical 
integrations ||. Two important questions remain, however, which the present section addresses. First, 
how badly does it fail when r c is not small enough to justify the Master-equation approach? Second, 
for realistic solar neutrinos what does the condition r c <C r osc translate into for the correlation length 
of a candidate fluctuation in the sun? 

Anwering the second of these first, for solar neutrinos the correlation time, r c , represents the 
time taken for a neutrino to traverse the correlation length, £, of a solar fluctuation. The numerical 
comparisons and analytical estimates of ref. |10[ indicate that the Master equation may be expected 
to work so long as £ satisfies £ < £ CT - lt = l/(eG F (n e )) (recall e 2 ~ (5n e d~n e ) / (n e ) 2 is a measure of 
the relative amplitude of the fluctuation at the position where the MSW resonance occurs). Taking 
e ~ 10% for the fluctuation amplitude and 4rit ~ 300km, corresponding to the electron density at the 
resonance point, this condition becomes £ < 3, 000km. 
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Figure 2: A comparison of the generalized Parke formula with the results of several numerical averages 
over a particular model - the 'Cell' model - of solar fluctuations. The survival probability is plotted as a 
function of the correlation length, I, using the representative values for solar neutrinos. See ref. for 
details. 

How badly does the Master equation describe fluctuations which are larger than this? Very badly 
indeed. Figure 2 shows a comparison of the generalized Parke formula with a series of direct numerical 
evaluations of neutrino evolution and the ensemble average, as reported in ref. |fL0|| . As may be seen 
from the figure, fluctuation effects completely disappear once £ passes appreciably above £ crit . 

There is a qualitative argument for why fluctuations with very long correlation lengths should 
not influence resonant neutrino conversion. A density fluctuation in the extreme limit £ ~ R Q , simply 



corresponds to an overall shift of the electron density throughout the entire sun. This has the sole 
effect of simply moving the position within the sun where the the resonance occurs, but does not alter 
the overall survival probability 

In summary, the lessons to be drawn from the extant calculations of neutrino oscillations within 
a fluctuating medium are these: 

1. Master equation methods describe solar neutrinos well for sufficiently small correlation lengths, 
but can dramatically overestimate the deviation from the MSW prediction when correlation 
lengths are large. For solar neutrinos, the dividing line between large and small correlation 
lengths is of order 3,000 km. 

2. Fluctuations can cause an appreciable deterioration of the MSW resonance only if two conditions 
are satisfied. First, fluctuation amplitudes cannot be too small. Second, their correlation lengths 
should be neither too big nor too small compared to the neutrino oscillation length. Most 
importantly, both of these conditions must apply at the positions in the sun where the neutrino 
passes through the MSW resonance. 

3. There is only one class of fluctuations in the sun which are known to exist with interestingly 
large amplitudes deep within the solar radiative zone, where the MSW resonance occurs. These 
are helioseismic waves |T2 ], for which the correlation lengths at resonance can be as large as 



Rq/10 ~ 70,000 km. As such they easily exceed the size of £ crit , and so preclude a description 
using only Master-equation techniques. 



3 Helioseismic Waves 

We now summarize the results of detailed simulations of neutrino oscillations in the presence of helio- 



seismic waves, as was originally reported in ref. |10|- This reference explicitly computes seismic wave 



profiles within the sun in order to reliably correlate their amplitude at the resonance point with ob- 
servations (which are made at the solar surface). It then numerically evolves neutrinos through these 
profiles to identify their effect on MSW oscillations. No use is made of the short-correlation-length 
approximation. 

There is a simple picture of what is going on. For any particular neutrino, microscopic fluctuations 
are irrelevant, and the perceived density profile is static since the time taken for a neutrino to traverse 
the sun (several seconds) is much shorter than the typical period of a helioseismic wave (minutes to 
hours). The oscillations of each neutrino are therefore well described by the ordinary Parke formula, 
i.e. eq. (|9|) with A = 1. The dominant change in the survival probability for each neutrino is then due 



to the differing scale heights, h, (and hence differing jump probabilities) which each neutrino sees. The 
average over the ensemble of density profiles seen by many neutrinos amounts to an average over h. 

These calculations find that helioseismic waves have no detectable effect on MSW oscillations, 
even under optimistic assumptions concerning their amplitudes. This conclusion turns out to rely in 
an important way on not using a Master-equation-based calculation, because those waves whose ampli- 
tudes can be large enough to affect the neutrino resonance, also have correlation lengths (wavelengths) 
which are much larger than £ CI .; t . More specifically: 

1. Those helioseismic waves (p- waves) which have already been observed at the solar surface gen- 
erally decrease in amplitude as one moves into the sun, with the result that their amplitude is 
too small in the MSW resonance region to produce observable effects. 

2. The main class of waves (g- waves) which are presently undetected are also too small to affect 
neutrino resonances if they carry energies which are similar to the p-waves which have been seen. 
This is so even though g-waves grow in amplitude with increasing solar depth. 

3. A particular class of g- waves is believed to be overstable, with runaway evolution once they 
are sufficiently excited. Although the endpoint of the resulting runaway behaviour remains 
controversial, it has been argued [1^] that the energy in these modes could be as much as 



10 10 times larger than what is in other modes. If the generalized Parke equation were used to 
propagate neutrinos through these waves, one would predict detectable deviations from MSW 
oscillations. Unfortunately, this prediction is incorrect since the wavelengths of the overstable 
modes are too large to permit a small-£ calculation. 
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